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摘 要: 针对 传统 狼 群 算法 (WPA) 存在 易 陷 入 局 部 最 优 解 、 计 算 资 源 耗 费 大 、 鲁 棒 性 低 等 问题 ， 提 出 一 种 基于 差分 
进化 的 改进 狼 群 算法 (DWPA)。 首 先 ， 通 过 引入 探 狼 搜索 因子 、 猛 狼 最 大 奔袭 次 数 、 自 适应 围攻 步 长 、 差 分 进化 策略 
等 对 传统 狼 群 算法 进行 了 改进 ， 降 低 算法 计算 耗费 的 同时 提高 了 算法 的 全 局 搜索 能 力 ; 然后 ， 运 用 马尔 可 夫 链 理论 证 
明了 DWPA 的 收敛 性 ; 最 后 ， 对 13 个 测试 函数 进行 寻 优 测试 并 与 WPA 等 4 种 算法 进行 对 比分 析 。 测 试 结果 表明 ， 
DWPA 具有 良好 的 彰 棒 性 和 全 局 搜索 能 力 ， 在 求解 多 峰 、 高 维 、 不 可 分 函数 方面 的 寻 优 能 力 尤为 突出 。 
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Research of improved wolf pack algorithm based on differential evolution 
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Abstract: Aiming at the problems of traditional wolf pack algorithm (WPA) : easy to fall in to local optimal, large computational 
resource cost and low robustness, propose an improved wolf pack algorithm based on differential evolution (DWPA) . First of 
all, propose search wolf search factor, maximum number of raid wolves, adaptive siege step size and differential evolution 
strategy to improve the traditional wolf pack algorithm, which can not only reduce the computational cost of the algorithm but 
also improve the global search ability. Then, prove the convergence of DWPA applying the Markov process. Finally, conduct 
optimization test on 13 functions and then compare it with WPA and other 4 algorithms. The test results show that DWPA has 
great robustness and global search ability, especially has an excellent optimizing ability in multi-peak, high-dimension, 
indivisible functions. 
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较 好 ， 已 被 广泛 应 用 于 无 人 机 航 迹 规划 中 、 水 电站 水 库 优 化 调 
度 03 等 人 类 生产 活动 中 。WPA 通过 模拟 狼 群 分 工 协作 捕猎 的 特 
群集 算法 主要 是 通过 模拟 生物 进化 和 生物 种 群 行为 来 求解 ” 征 ， 抽 象 出 游 走 、 召 唤 、 奔 袭 和 围攻 等 行为 与 “ 胜 者 为 王 ” 的 


优化 问题 的 智能 算法 卜 。 常 见 的 群集 智能 算法 主要 有 粒子 群 算 头 狼 产 生机 制 和 “ 强 者 生存 ”的 群体 更 新 机 制 来 进行 优化 问题 


(PSO) 2-4、 和 鱼 群 算法 (FSA) G5.5]、 蚁 群 算法 (ACO) 中、 求解 ， 具 有 较 好 的 寻 优 性 能 。 但 也 不 可 避免 地 存在 一 些 不 足 之 


差分 进化 算法 (DE)9 和 等， 它们 在 解决 优化 问题 时 各 有 千秋， 为。 处 易 陷入 局 部 最 优 解 、 计 算 资 源 耗 费 大 、 和 鲁 棒 性 低 等 。 文 册 
某 些 复杂 问题 的 解决 提供 了 良好 的 思路 。 [13] 提 出 一 种 改进 的 狼 群 算法 ， 该 算法 根据 传统 狼 群 算法 基本 
狼 群 算法 最 早 是 由 Yang 等 人 外 提出 ，2013 年 ， 吴 虎 胜 等 ”思想 提出 探 狼 更 新 规则 并 引入 相位 因子 ， 同 时 优化 了 传统 狼 群 
和 人 09 在 分 析 狼 群 分 工 协作 围捕 猪 物 行 为 的 基础 上 ， 提 出 了 一 种 算法 步 长 的 种 类 ， 设 计 了 新 的 猛 狼 位 置 更 新 公式 ， 通 过 测试 函 


\ 


全 新 的 狼 群 算法 (WPA)。 算法 虽 问 世 时 间 较 短 , 但 由 于 其 性 能 。 数 仿 真 模 拟 验 证 了 该 算法 的 有 效 性 。 文 献 [14] 结 合 文化 算法 提 
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录用 稿 


出 了 一 种 文化 狼 群 算法 ， 该 算法 可 有 效 解决 人 工 狼 搜索 的 盲目 
性 问题 ， 通 过 对 三 个 复杂 函数 的 测试 分 析 ， 验 证 了 该 文化 狼 群 
算法 的 有 效 性 。 文 献 [15] 将 PSO 算法 中 求解 当前 局 部 最 优 的 思 
想 引 入 到 狼 群 算法 的 游 走 和 召唤 行为 中 ， 并 利用 混沌 法 对 得 到 


rand e(0,1) 内 的 随机 数 。 


2) 头 狼 产生 规则 
在 初始 化 产生 的 种 群 中 ,通过 式 (2) 计算 每 一 区 人 工 狼 的 


的 次 优 解 进行 优化 。 改 进 后 的 算法 大 大 提高 了 搜索 的 准确 度 
避免 了 陷入 局 部 极 值 ， 并 通过 仿真 验证 了 该 算法 的 有 效 性 。 


适应 度 值 ， 适 应 度 值 越 小 的 位 置 视 为 猎物 浓度 越 大 ， 距 离 所 求 
最 优 解 越 近 ， 选 取 有 具有 最 优 适 应 度 值 的 人 工 狼 作 为 本 次 迭代 过 


针对 WPA 存在 的 不 足 之 处 ， 本 文 提出 一 种 基于 差分 进化 
的 改进 狼 群 算法 (DWPA )。 该 算法 在 对 传统 WPA 算法 智能 行 
为 改进 的 基础 上 , 引入 差分 进化 中 的 交叉 、 变 异 、 选 择 等 操作 ， 
大 大 提高 了 算法 的 寻 优 性 能 及 鲁 棒 性 。 


1 ”DWPA 算法 原理 


DWPA 算法 继承 了 WPA 算法 的 所 有 行为 并 对 其 进行 了 改 
进 。 主 要 改进 如 下 : a)WPA 中 游 走行 为 缺乏 引导 ， 一 旦 游 走 方 
向 数 h 确定 ， 游 走 方向 便 随 之 确定 ， 无 法 全 面 搜索 解 空间 ， 易 
使 算法 陷入 局 部 最 优 , 故 在 DWPA 游 走行 为 中 引入 探 狼 搜索 因 
子 ,使 得 算法 寻 优 能 遍历 整个 空间 ;b)WPA 中 的 奔袭 行为 需要 分 
析 计 算 人 工 狼 与 头 狼 之 间 的 距离 并 要 求 所 有 人 工 狼 与 头 狼 之 间 
距离 小 于 临界 距离 时 才 结 束 奔袭 行为 ， 在 面 对 高 维 函数 的 寻 优 
时 ， 必 然 造成 大 量 计算 资源 的 耗费 ， 同 时 由 于 临界 距离 这 一 唯 
一 行为 终止 条 件 的 存在 将 可 能 使 得 算法 无 法 跳出 奔 缆 行为 ， 降 
低 算 法 的 鲁 棒 性 , 故 在 DWPA 奔袭 行为 中 除去 人 工 狼 距离 计算 
过 程 ， 并 引入 猛 狼 最 大 奔袭 次 数 和 头 狼 蔡 换 两 个 条 件 作 为 行为 
终止 条 件 ， 防 止 算法 陷入 无 限 循环 并 提高 计算 效率 和 算法 鲁 棒 
性 ;c)WPA 中 国 攻 行为 不 具有 自 适应 特性 , 决策 变量 取 值 范围 和 
步 长 因子 一 旦 确定 ， 所 有 人 工 狼 将 以 同一 确定 步 长 进行 围攻 ， 
故 在 DWPA 围攻 行为 中 引入 自 适 应 围攻 步 长 , 使 得 距离 头 狼 相 
对 较 远 的 人 工 狼 以 相对 较 大 的 步 长 向 头 狼 靠 近 ， 距 离 头 狼 相对 
较 近 的 人 工 狼 以 较 小 步 长 对 头 狼 附 近 区 域 进行 更 加 精细 的 搜索 ， 
兼顾 了 算法 的 局 部 搜索 和 全 局 搜索 能 力 。 同 时 ， 针 对 WPA 算 
法 本 身 易 陷 入 局 部 最 优 的 问题 ， 结 合 差分 进化 算法 在 局 部 寻 优 
方面 的 优势 外, 在 改进 WPA 算法 中 加 入 差分 进化 的 变异 、 交 又 
和 选择 操作 ,使 算法 能 跳出 局 部 最 优 ,提高 算法 全 局 搜索 能 

综 上 ，DWPA 算法 由 游 走 、 召 唤 、 奔 袭 、 围 攻 、 变 异 、 交 
又 、 选 择 等 七 种 行为 及 “ 胜 者 为 王 ” 的 头 狼 产 生 、“ 强 者 生存 ” 
的 群体 更 新 等 两 种 机 制 组 成 。 

下 面 以 最 小 值 寻 优 问题 为 例 , 详细 
过 程 : 

1) 狼 群 初始 化 

采用 式 (1) 进行 狼 群 初始 化 ， 产 生 n 匹 人 工 狼 : 


上 


介绍 DWPA 算法 的 主要 


XI = 和 -+ rand x (XY 一 =1,2,…,n (1) 


其 中 : X; 表示 第 t 代 种 群 中 第 i 还 人工 狼 ; X* 和 XX" 分别 为 变 


邱 


=(w,b,…,xp) 取 值 的 下 界 和 上 界 ，D 为 变量 维 数 ， 


程 中 的 头 狼 ， 若 有 多 匹 人 工 狼 适 应 度 值 都 是 最 优 值 ， 则 任 选 一 
匹 作 为 头 狼 ， 头 狼 适 应 度 值 为 2 ; 在 进化 过 程 中 ， 将 每 次 进 
化 后 具有 最 优 适应 度 值 的 人 工 狼 与 前 一 代 中 的 头 狼 进行 比较 ， 
若 当前 人 工 狼 的 适应 度 值 更 优 ， 则 该 人 工 狼 蔡 代 头 狼 成 为 新 的 
头 狼 。 对 于 最 小 值 问 题 ， 当 QQ < Qi 时 ， 视 为 8 优 于 Qon 。 


Q=f(X) 0) 


其 中 : f(X) 为 所 求 最 小 值 问 题 的 目标 函数 。 
3) 游 走行 为 
除 头 狼 外 ， 在 当前 猎物 群 中 选取 适应 度 值 最 优 的 Sm 匹 人 
工 狼 作 为 探 狼 ， 在 空间 中 搜寻 猎物 。 sm 取 值 为 


[ay(c+D,n/o] 间 的 整数 , 其 中 ，c 为 探 狼 比例 因子 。 探 狼 朝 


向 h 个 方向 分 别 前 进一步 并 记录 每 前 进一步 后 所 处 位 置 的 适应 
度 值 ， 然 后 退回 原 位 置 ， 则 向 第 pp=1,2,…,D) 个 方向 前 进 后 ， 探 
狼 i 在 第 d(d=1,2,…,D) 维 空间 中 所 处 的 位 置 如 式 〈3 ) 所 示 。 
x = Xiy + yx step’ (3) 

其 中 : ”step, 为 探 狼 游 走 步 长 ，7 为 探 狼 搜索 因子 ， 为 [-1,1 内 
的 随机 数 。 

此 时 , 探 狼 所 处 位 置 的 适应 度 值 为 Qip, 选择 适应 度 值 最 小 
且 小 于 当前 位 置 适 应 度 值 Qio 的 方向 前 进一步 ， 更 新 探 狼 信息 
Xi。 将 一 次 游 走行 为 结束 后 的 具有 最 小 适应 度 值 Qmin 的 人 工 狼 
与 头 狼 Qiead 进行 比较 ， 若 Qmin<Qiead， 则 Qiead=Qmnin， 适 应 度 
最 小 的 探 狼 蔡 代 原 头 狼 成 为 新 的 头 狼 , 并 发 起 召唤 行为 ; 否则 ， 
重复 游 走行 为 ,， 直到 某 匹 狼 i 所 处 的 位 置 适应 度 值 Q<Qead， 
游 走 次 数 Ti 达到 最 大 游 走 次 数 Timax。 需 要 指出 的 是 ， 每 匹 探 
狼 由 于 适应 度 值 存在 差异 ， 其 搜索 方式 也 存在 差异 ， 故 每 匹 探 
狼 h 的 取 值 是 不 同 的 ， 将 在 [hmin,hmax] 内 取 随 机 整数 ， 一 般 来 说 
h 越 大 ， 探 狼 搜 索 越 精细 但 也 将 影响 算法 运行 速率 。 

4) 奔袭 行为 

头 狼 通 过 喷 叫 发 起 召唤 行为 ， 召 集 附近 狐 狼 快速 向 头 狼 靠 
近 。 猛 狼 i 在 第 kt+l 次 进化 时 ， 在 第 d(d=1,2,…,D) 维 空间 中 所 
处 的 位 置 为 如 式 〈4) 所 示 。 


mr 


大 十 1 
Xa 二 


k a kk kk 
Xia + Step, x(g’ - 蕊 )/|g’ Xa 


(4) 


其 中 : step， 为 独 狼 奔袭 步 长 ， g4 为 第 k 代 群 体 头 狼 在 第 
d(d=1,2,…,D) 维 空间 中 的 位 置 。 

在 狐 狼 奔袭 的 过 程 中 ， 将 每 次 奔袭 结束 后 县 有 最 小 适应 度 
值 Quan 的 猛 狼 与 头 狼 Qiead 进 行 比 较 , 若 Qmin<Qiead; 则 Qiead=Qmin， 
该 狐 狼 成 为 新 的 头 狼 ， 奔 袭 行为 结束 ， 进 入 下 一 行为 ， 否 则 继 
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续 奔袭 直到 某 匹 猛 狼 i 的 适应 度 值 小 于 头 狼 适应 度 值 或 奔袭 达 
到 最 大 奔袭 次 数 Tzmax。 

5) 围攻 行为 

猛 狼 奔袭 行为 结束 后 ， 将 联合 所 有 的 探 狼 ， 以 头 狼 所 在 的 


CR 越 小 越 有 利于 保持 狼 群 多 样 性 和 全 局 搜索 。 
8) 选择 行为 
狼 群 经 过 交叉 行为 后 ， 对 于 得 到 的 人 工 狼 试验 矢量 

必 ” 与 人 工 狼 目标 矢量 Xi , 若 试验 矢量 的 适应 度 值 小 于 目标 矢 


位 置 为 猎物 移动 位 置 ， 执 行 围 攻 行 为 。 对 于 第 k 代 狼 群 ， 假 定 
头 狼 在 第 d(d=1,2,.…,D) 维 空间 中 的 位 置 为 G6;, 则 狼 群 的 围攻 行 
为 如 式 (5) 所 示 。 


Xi = x +Axstep" x(G* 一 %)/|G’ 一 项 | (5) 


其 中 : 4 为 [-1,1] 间 均匀 分 布 的 随机 数 ，step; 为 人 工 狼 i 执 行 盏 
攻 行 为 时 的 自 适 应 围攻 步 长 ， 距 离 头 狼 相对 较 远 的 人 工 狼 以 相 
对 较 大 的 步 长 向 头 狼 靠近 ， 距 离 头 狼 相 对 较 近 的 人 工 狼 以 较 小 
步 长 对 头 狼 附 近 区 域 进 行 更 加 精细 的 搜索 。 以 搜索 空间 上 下 界 
和 搜索 空间 中 的 适应 度 最 差 值 与 最 优 值 的 差 作为 基准 值 ， 以 自 
适应 围攻 步 长 与 人 工 狼 适 应 度 值 和 头 狼 适应 度 值 的 差 成 正 相 关 
关系 来 建立 自 适 应 围攻 步 长 step; 如 式 〈6) 所 示 。 
step,= (OQ,-0,, Dr (6) 

其 中 : max 0, 和 min @ 分 别 为 算法 初始 化 产生 的 人 


狼 的 最 差 


量 , 则 用 wu 取代 Xt* 成 为 第 k+1 代 并 更 新 人 工 狼 位 置 , 否则 ， 
直接 用 和 作为 第 krl 代 更 新 ， 即 式 (10)。 
人 
Xk Ql(u'") > Q(x) 


对 于 得 到 的 人 工 狼 矢量 X*”， 找 出 其 中 具有 最 小 适应 度 值 
Qnmin 的 人 工 狼 并 与 头 狼 适 应 度 Qiead 进行 比较 ， 若 Qmin<Qiead， 
则 Quead=Qmin， 适 应 度 值 最 小 的 人 工 狼 选 为 头 狼 。 

9) 狼 群 更 新 机 制 

狼 群 捕食 得 到 的 猎物 ， 并 不 是 平均 分 配给 每 一 匹 狼 ， 将 按 
照 “ 肉 弱 强 食 ” 的 原则 进行 食物 分 配 ， 这 就 可 能 导致 弱小 的 狼 
被 饿 死 ， 这 一 原则 通过 在 算法 中 去 掉 适 应 度 值 最 差 的 R 匹 狼 ， 
同时 按 式 (1) 产生 R 匹 人 工 狼 来 进行 模拟 。R 为 


[n/(2xB),n1B] 之 间 的 随机 整数 ，B 为 狼 群 更 新 比例 因子 。 


通过 以 上 行为 ， 使 得 狼 群 不 断 进 行 更 新 进化 ， 直 到 算法 达到 限 


适应 度 值 和 最 优 适 应 度 值 。 
若 围 攻 后 的 人 工 狼 所 处 位 置 的 适应 度 值 小 于 原 位 置 的 适应 
度 值 ， 则 更 新 人 工 狼 位 置信 息 ， 否 则 ， 人 工 狼 位 置 不 变 。 
特别 地 : 游 走 步 长 stepa, 奔袭 步 长 stepe 在 待 寻 优 空间 中 存 
在 如 式 〈7) 所 示 关 系 : 


step, = step, 12=|X —X"|/S (7) 


其 中 : S 为 步 长 因子 ， 表 示人 工 狼 在 待 寻 优 空间 中 搜索 的 精细 
程度 。 

6) 变异 行为 
在 执行 围攻 行为 后 的 群体 中 ， 随 机 选择 两 匹 不 同 的 人 工 狼 
(第 k 代 ) Xr 、X! ,将 两 者 的 矢量 差 赋予 权 值 后 与 第 三 匹 随 
机 选择 的 人 工 狼 X* 相 加 , 得 到 人 工 狼 变异 矢量 vw ,对 于 任意 
人 工 狼 目 标 矢量 X; 有 


WW = XI+F* (XX) (8) 


其 中 : F 为 缩放 因子 ， 一般 在 [0,2] 中 取 值 ， 表 示 对 矢量 差 的 缩 
放 程 度 。 

7) 交叉 行为 

狼 群 经 过 变异 行为 后 ， 对 于 得 到 的 人 工 狼 变异 矢量 站 ”和 
人 工 狼 目标 矢量 XX  ， 按 式 (9) 进行 参数 混合 ， 生 成 人 工 狼 试 


验 矢 量 必 ” : 


x rand(j) > CRorj # randn(i) 


el | rand(}) < CRorj = randn(i) , 
y 
其 中 : j 表示 人 工 狼 的 第 j 个 变量 ，randn(i) e[l,2,…,D]，, 为 随 
机 选择 的 维 数 变量 的 标识 ;CR 为 交叉 概率 因子 , 一 般 在 [0,1] 中 


取 值 ，CR 越 大 越 有 利于 算法 进行 局 部 搜索 和 加 快 收敛 速度 ， 


定 的 最 大 迭代 次 数 G 或 算法 获得 的 最 优 值 达 到 预 设 精度 时 结 
束 ， 故 有 DWPA 求解 优化 问题 流程 图 如 图 1 所 示 。 


初始 化 狼 群 位 置 Xi 


计算 各 人 工 狼 适 应 
度 值 Qi 


| 


一 


根据 头 狼 产生 规则 确定 maxQo 和 Ce 
确定 头 狼 min 狼 群 执行 变异 行为 
探 狼 执行 游 走行 为 狼 群 执行 交叉 行为 


狼 群 执行 选择 行为 


更 新 当前 头 狠 位 Ce 

置 ， 并 发 起 召唤 要 新 次 急 位 攻 
| 执行 “ 强 者 生存 ” 

猛攻 执行 友 交 行为 ye 


输出 头 狼 位 置 ， 即 
为 所 求 最 优 解 


更 新 头 狼 位 置 


狼 群 执行 围攻 行为 


1 DWPA 求解 优化 问题 流程 图 


录用 稿 


2 DWPA 收敛 性 分 析 


马尔 可 夫 (Markov) 链 是 一 种 无 后 效 性 的 随机 过 程 ， 常 被 


用 于 i 
歼 、 


FE 明 算法 的 收敛 性 
围攻 、 交 义 、 变 异 、 
每 种 行为 都 只 与 当前 群体 状态 有 关 而 与 之 前 的 状态 
可 知 ，DWPA 种 群 序列 符合 Markov 链 。 


0194。DWPA 是 不 断 重 复 游 走 、 召 唤 、 奔 
选择 和 狼 群 更 新 行为 的 过 程 ， 种 群 的 
无 关 ， 由 此 


引 理 1 文献 [17] 
件 :a) 对 解 空 
子 可 达 的 ; b) 种 群 


序列 


zs 间 中 任意 两 解 丸和 


若菜 智 


下 
2X2 古 略 


明 j 能 算法 满足 如 下 两 个 条 
算法 中 的 各 种 算 


0O2COv 是 单调 的 。 则 该 智能 算法 以 


中 证 


概率 1 收敛 于 问题 的 全 


局 最 优 解 。 


通过 以 上 结论 可 知 , 如 果 DWPA 能 够 满足 引 理 1 的 两 个 条 
件 ， 则 说 明 DWPA 是 以 概率 1 收敛 问题 的 全 局 最 优 解 的 。 
参考 引 理 1 的 结论 , 对 于 DWPA 收敛 性 问题 , 若 其 满足 以 
下 两 个 条 件 :a)DWPA 种 群 序列 的 Markov 链 是 遍历 的 ;b)DWPA 
吕 的 序列 解 是 一 个 有 限 齐 次 Markov 链 。 则 DWPA 将 以 概率 1 收 
~ 敛 于 问题 的 全 局 最 优 解 。 
© 假设 DWPA 搜索 的 解 空间 为 Q , 游 走 、 奔 袭 、 围 攻 、 交叉 、 
EC 变异 、 选 择 和 狼 群 更 新 行为 所 引起 的 状态 空间 的 状态 转移 分 别 
加 用 转移 矩阵 人 B、 从 小 从 XX 和 6 表示 , 则 DWPA 种 群 的 Markov 
< 十 链 的 转移 矩阵 定义 为 
©O P=YxBxWxYJxVxXxG (1) 
a 在 证 明 DWPA 的 收敛 性 之 前 ， 首 先 给 出 以 下 几 个 定义 : 
ew _ 
EN 定义 1 假设 为 Markov 链 的 转移 概率 矩阵 ， 若 有 对 于 
7 vvjso, 尖 >1, 使 得 五 >0, 则 称 该 Markov 链 是 不 可 约 的 。 
全 定义 2 假设 有 非 空 集合 1={K|k>Z 太 >0,Vi,jeQ}， 若 
pT 该 集合 7 的 最 大 公约 数 为 1 ， 则 称 该 Markov 链 是 非 周期 的 。 
定义 3 假设 所 = 并 ,对 于 常 返 状态 记 若 有 从 <+oo， 
则 称 & 为 正常 返 的 。 特别 地 , 当 i 为 正常 返 的 且 非 周期 的 , 则 该 
Markov 链 是 遍历 的 。 
押 证 明 DWPA 收敛 性 : 
推理 1 DWPA 种 群 序列 的 Markov 链 是 遍历 的 。 
推理 1 的 证 明 ， 将 分 以 下 三 个 步骤 进行 
引证 明 DWPA 的 种 群 序列 Markov 链 是 不 可 约 的 。 


假设 Q. = {X, xX 


Xi; 为 第 i 


X、} 为 算法 的 第 k 代 种 群 ， 


匹 人 工 狼 的 状态 。 


于 种 群 Markov 链 的 转移 概率 矩阵 


妃 =P{O= 


Q. >0 恒 成 立 ， 则 有 种 群 
序列 的 Markov 链 是 不 可 约 的 。 


1，DWPA 种 群 


中 Q. =i,k>]} 内 与 始末 状态 ij 相关 ， 


转移 概率 矩阵 P 为 


E 定 矩阵 ， 由 定义 


b) 证 


明 DWPA 种 群 


的 Markov 链 是 非 周期 且 


不 可 允 


王 盈 祥 


于 某 一 给 定 的 k>0， 
可知 ， 必 然 忆 Ee 0， 


链 是 非 
且 不 可 约 
C) 订 


期 的 ， 
的 。 


,为 状态 i 经 过 各 种 行为 转化 为 状态 j 的 概率 ， 


， 等 : 基于 差分 进 


由 


hi Nn aXI 上 


进化 的 改选 狠 


日 


明和 


法 研究 


(不 


二 


得 到 的 算法 Markov 链 是 不 


E 明 DWPA 种 群 的 Markov 链 是 遍历 的 。 
由 于 转移 矩阵 YB， 人 办 人 内 XX 和 6 的 值 均 在 [0,1 内， 且 


使 得 驴 >0 成 立 ， 结 合 定义 2 可 知 
=1 。 因 此 可 得 I 的 最 大 公约 数 为 1, 故 有 DWPA 种 群 
周 结合 四 可 知 DWPA 种 群 的 Markov 链 是 非 周 期 


的 Markov 


必 有 0<P<1， 


设 o=max{P, [vi,j eQ} ， 则 由 Canchy-Riemann 方程 及 定义 3 


可 知 : 4 = DP sy kat<a。 
k k 


综 上 分 析 可 知 ，DWPA 种 群 序列 上 


理 1 得 证 。 


推理 2 DWPA 得 序列 
证 明 


解 


又 、 变 蜡 、 
个 体 的 更 新 都 具有 优质 选择 的 


于 每 一 代 种 群 Q. 都 是 有 
是 有 限 的 ， 其 次 ， 算 法 不 断 重复 游 走 、 召 唤 、 奔 袭 、 围 攻 、 
选择 和 狼 群 更 新 行为 的 过 程 寻找 更 加 优质 的 猎物 ， 

群体 第 kt+1 代 的 产生 仅 与 


个 有 限 齐 次 


了 


特性 。 


限 的 ， 故 


的 Markov 链 是 遍历 的 ， 


Markov 链 
其 Markorv 链 必 


RE 


父 


第 k 代 有 关 , 而 与 ,和 kk 无关 。DWPA 经 过 反复 迭代 寻 优 ， 必 


能 得 到 一 组 序列 解 ， 且 该 序列 
里 2 得 证 。 


上 推理 


综 


个 有 限 齐 次 


Ee 


以 概率 1 收敛 于 所 求 问 题 的 全 局 最 优 解 。 


3 ”DWPA 算法 仿真 测试 


3.1 基本 测试 函数 介绍 


为 了 测试 DWPA 算法 的 怕 


分 析 


能 及 其 有 效 性 ， 


单 


数学 特征 和 规模 不 同 的 标准 
测试 函数 如 表 1 所 示 ; 在 表 1 


标 测 试 函数 对 
征 栏 中 


特 和 


Markov 链 ， 推 


1 和 2 可 知 ，DWPA 种 群 序 列 的 Markov 链 是 遍 
历 的 且 其 序列 解 是 一 个 有 限 齐 次 Markov 链 ， 至 此 证 


明 DWPA 


本 文选 用 13 个 
其 进行 测试 ， 各 


，“U” 表 示 单 峰 函 数 ， 


即 该 函数 在 定义 域内 只 有 一 个 极 值 ， 即 全 局 最 优 


最 优 值 ;“M” 表 示 多 峰 函 
极 值 , 主要 用 来 检测 算法 


的 全 局 


值 ， 没 有 局 部 


即 该 函数 在 定义 


寻 优 和 避免 早熟 的 能 力 08， 


域内 有 多 个 局 部 
“Gg” 


表示 可 分 函数 ,，“N” 
数 能 ) 


表 不 个 可 


分 函数 ， 若 某 具 有 mn 


个 变量 的 函 


j 这 mn 个 单 变量 函数 之 和 来 表示 该 函数 ， 则 该 函数 为 可 分 


函数 ， 否 则 ， 为 不 可 分 了 
复杂 ， 在 寻 优 过 程 中 对 算法 的 
栏 中 ， 集 合 部 分 为 变量 取 值 范 


数 。 不 可 分 函数 相对 本 


性 能 要 求 更 高 。 
围 ， 上 标 为 函 


数 


也 越 高 。 
3.2 算法 对 比 验证 及 分 析 
为 了 测试 DWPA 算法 的 性 


能 ,本 文 同 时 用 


分 函数 将 更 加 
在 表 1 的 定义 域 
的 维 数 D， 测 试 


售 数 DD 越 高 ， 算 法 的 运算 复杂 度 越 大 ， 对 算法 性 能 的 要 求 


DWPA、WPA、 


:201804.02401v1 


chinaXiv 
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PSO、DE 和 FSA 算法 对 以 上 13 个 函数 分 别 重复 进行 100 次 寻 


4590 CPU @3.30GHz， 内 存 16 GB，Win7 专业 版 64 位 操作 系 


优 测试 并 记录 寻 优 的 最 优 值 、 最 差 值 、 平 均值 、 标 准 差 和 计算 统 ,，MATLAB R2013b。 
成 功率 五 个 方面 进行 对 比分 析 。 实 验 环境 为 : Intel® Core™Mi5- 
表 1 测试 函数 
、 ps 上 局 
函数 名 函数 表达 式 特征 定义 域 
最 优 
2 学 
Sphere A US [-10,10]? 0 
i=] 
Matyas f=026(® + )-0.48n% UN [-10,10] 0 
Booth ={w t= MS [-10,10]? 0 
Eggcrate f= 其 + 并 + 25(sin™x, + Sin?x, ) MN [-x, xz] 0 
D D 
Trid6 k= 2 .00-17 -Ys UN [-36,36]s -50 
让 了 这 2 
区 党 
Sumsquares-1 f= US [-10,10]" 0 
这 了 
区 i 
Rastrigin 万 => Lx -10cos(2zx; )+10] MS [-10,10]60 0 
这 了 
Griewank A= $e Ie) + MN 600,600]1% 0 
riewanl 7 000 和 ”| 所 [-600,600] 
D 和 
Quadric 万 => (2 xD) MS [-30,30]1% 0 
i=1 k=l 
页 
Sumsquares-2 f= US [-10,10]'" 0 
j=] 
oe 和 记 
Rosenbrock f= [100 -BP +-17] MN [-30,30]1% 0 
i=] 
1 ， 7 了 忆 
Ackley-1 fi = -20exp(-0.2 万 之 二 )- ep{ F202 ))+20+e MN [-32，32]100 0 
el 让 
7 忆 1 己 
Ackley-2 fi = -20exp(-0.2 万 之 刀 )-ep( F202 ))+20+e MN [-32，32]200 0 
i=J el 


计算 成 功率 定义 为 成 功 计 算 次 数 占 总 计算 次 数 的 百分比 
09。 对 于 每 次 寻 优 计 算得 到 的 结果 Q， 当 Q 满足 式 (12) 时 ， 
表示 此 次 寻 优 成 功 。 


0 CO 
best 


C 一 CO， 


其 中 : 


Qs 为 所 求 函数 的 到 


<1E-03 O zx0 


<1E-03 Q,,=0 


表 2 各 算法 参数 对 照 表 


(12) 


E 论 全 局 最 优 值 。 


寻 优 计算 的 最 大 迭代 次 数 设 为 2000, DWPA、WPA、PSO、 
DE 和 FSA 算法 的 种 群 规模 设 为 50， 其 中 : WPA 算法 的 参数 


依据 文献 [10] 选 取 , PSO 算法 的 参数 依据 文献 [2] 选 取 , DE 算法 
的 参数 依据 文献 [8] 选 取 , FSA 算法 的 参数 依据 文献 [6] 选 取 , 整 


里 如 表 2 所 示 ; 五 种 算法 对 以 上 13 个 函数 寻 优 计算 的 结果 如 


表 3 所 示 。 


算法 参数 取 值 
惯性 权重 w=0.7298 ， 学 习 因 子 c = c =2 ， 个 体 速 度 范 
PSO 
围 [-0.5,0.5] 
DE 缩放 因子 F-0.6， 交 叉 概 率 因子 CR=0.9 
最 多 试探 次 数 trynum=100, 感 知 距离 =1， 拥 挤 度 因子 
人 de=0.618， 步 长 step=0.1 
探 狼 比例 因子 <=4 ， 最 大 游 走 次 数 Taax=20, 距 离 判定 因 
子 w=500 ， 步 长 因子 s=1000, 更 新 比例 因子 B=6 
探 狼 比例 因子 a=4 ， 最 大 游 走 次 数 Timax=10, 最 大 奔袭 次 
DWPA ” 数 Tna=10, 步 长 因子 S=1500, 更 新 比例 因子 p=3， 缩 放 医 


子 F=0.6， 交 又 概 率 


因子 CR=0.9 


Matyas,5 种 算法 都 寻 优 成 功 且 具 


a) 对 于 单 峰 、 


表 3 中 各 种 算法 的 寻 优 计算 结果 可 以 看 出 : 


低 维 的 可 分 函数 Sphere 及 不 可 分 函数 


有 较 好 的 改 


FE 能 , 其 中 DWPA 和 


(I | 
dd 作 
录用 稿 王 盈 祥 ， 等 : 基于 差分 人 合作 部 


WPA 的 寻 优 能 力 相当 ， 尤 其 在 求解 不 可 分 函数 Matyas 时 ， 功 且 性 能 较 好 , DE 和 WPA 存在 陷入 局 部 最 优 而 寻 优 失 败 的 情 
DWPA 和 WPA 性 能 远 高 于 其 他 三 种 算法 。 况 ，FSA 寻 优 完全 失效 。 
b) 对 于 多 峰 、 低 维 的 可 分 函数 Booth，5 种 算法 都 寻 优 成 功 e) 对 于 单 峰 、 可 分 函数 Sumsquares-1, DWPA、WPA、PSO、 
且 性 能 都 较 好 ， 其 中 PSO 寻 优 性 能 最 佳 ， 此 时 DWPA 性 能 优 。 FSA 都 寻 优 成 功 , 其 中 DWPA 和 WPA 具有 最 佳 寻 优 性 能 , D 
势 并 不 明显 。 寻 优 完全 失效 。 对 于 高 维 函 数 ，PSO、DE 和 FSA 都 无 法 跳出 
c] 对 于 多 峰 、 低 维 的 不 可 分 函数 Eggcrate,FSA 存在 由 于 陷 ”局 部 最 优 而 导致 寻 优 失败 ， 不 能 获得 较 好 的 寻 优 效果 只 有 


入 局 部 最 优 值 而 使 寻 优 失败 的 情况 ， 使 其 寻 优 成 功率 较 低 ， 其 “WPA 和 DWPA 有 一 定 的 寻 优 能 力 ， 故 以 下 只 对 比 WPA 和 
他 四 种 算法 全 部 寻 优 成 功 ， 此 时 DWPA 具有 最 佳 寻 优 性 能 。 DWPA 的 寻 优 性 能 。 
d) 对 于 单 峰 、 不 可 分 函数 Tiid6，DWPA 和 PSO 都 寻 优 成 
表 3 各 算法 寻 优 结果 


函数 算法 最 优 值 最 差 值 平均 值 标准 差 计算 成 功率 
PSO 5.87E-108 2.51E-99 2.35E-100 6.30E-100 100% 
DE 2.37E-23 6.24E-16 5.88E-18 8.38E-16 100% 
Sphere 
FSA 2.40E-16 4.73E-09 8.64E-10 7.53E-09 100% 
(2 维 ) 
WPA 0 0 0 0 100% 
DWPA 0 0 0 0 100% 
PSO 2.33E-13 9.70E-10 1.11E-11 2.01E-11 100% 
DE 1.97E-17 9.77E-15 5.01E-16 2.31E-16 100% 
Matyas 
FSA 8.33E-14 5.78E-11 9.89E-12 1.32E-11 100% 
(2 维 ) 
WPA 3.01E-101 6.18E-94 2.97E-96 1.12E-94 100% 
DWPA 2.23E-101 5.65E-95 3.42E-97 7.20E-96 100% 
PSO 0 0 0 0 100% 
DE 6.55E-20 1.87E-17 2.56E-18 8.11E-17 100% 
Booth 
FSA 1.67E-11 2.43E-09 5.98E-10 6.27E-10 100% 
(2 维 ) 
WPA 7.22E-10 6.46E-07 2.13E-08 1.71E-08 100% 
DWPA 6.15E-11 5.45E-09 1.12E-10 1.33E-09 100% 
PSO 2.52E-189 9.47E-182 5.21E-183 8.57E-183 100% 
DE 4.31E-22 6.66E-19 1.85E-19 1.74E-19 100% 
Eggcrate 
FSA 5.74E-07 9.49E+00 1.52E+00 3.51E+00 5% 
(2 维 ) 
WPA 1.27E-171 1.25E-164 2.65E-165 0 100% 
DWPA 0 0 0 0 100% 
PSO -50 -50 -50 0 100% 
DE -50 -49.9387 -49.9969 0.0134 97% 
Trid6 
FSA -49.9498 -49.9395 -49.9476 0.0041 0 
(6 维 ) 
WPA -50 -49.9497 -49.9975 0.011 96% 
DWPA -50 -50 -50 0 100% 
PSO 2.21E-46 1.16E-34 5.92E-36 2.59E-35 100% 
DE 2.53E-02 1.34E-01 7.58E-02 2.11E-02 0 
Sumsquares-1 
FSA 2.13E-12 1.23E-07 2.12E-09 1.87E-09 100% 
(10 维 ) 
WPA 0 0 0 0 100% 
DWPA 0 0 0 0 100% 
PSO 1.80E+02 2.68E+02 2.25E+02 2.18E+01 0 
DE 4.80E+02 5.78E+02 5.33E+02 1.99E+01 0 
Rastrigin 
FSA 1.10E+02 1.22E+03 9.74E+02 5.85E+01 0 
(60 维 ) 
WPA 4.41E-09 3.37E-06 3.75E-07 7.48E-07 100% 


DWPA 0 2.41E-12 1.58E-13 5.20E-13 100% 
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PSO 1.79E-01 2.55E-01 2.21E-01 2.06E-02 0 
DE 3.27E-02 9.19E-01 1.01E-01 2.14E-02 0 
Griewank 
FSA 7.83E+03 1.01E+04 9.83E+03 4.54E+02 0 
(100 维 ) 
WPA 3.11E-12 3.88E-01 1.49E-01 9.73E-02 94% 
DWPA 2.44E-15 6.68E-12 1.59E-12 1.95E-12 100% 
PSO 4.37E+02 1.13E+03 7.04E+02 1.50E+02 0 
DE 1.19E+05 2.30E+05 1.83E+05 3.29E+04 0 
Quadric 
FSA 1.31E+05 3.50E+05 2.78E+05 5.56E+04 0 
(100 维 ) 
WPA 5.93E-10 8.64E-07 2.33E-07 3.06E-07 100% 
DWPA 4.15E-12 8.48E-08 1.11E-08 1.89E-08 100% 
PSO 5.98E+02 1.14E+03 8.57E+02 1.34E+02 0 
DE 5.73E-01 4.72E+03 2.61E+03 1.23E+03 0 
Sumsquares-2 
FSA 3.03E+05 4.16E+05 3.61E+05 2.39E+04 0 
(100 维 ) 
WPA 5.02E-09 5.02E-07 1.21E-07 1.37E-07 100% 
DWPA 7.15E-16 6.12E-12 8.78E-13 1.43E-12 100% 
PSO 1.16E+01 1.45E+01 1.29E+01 5.83E-01 0 
DE 5.11E+00 1.58E+01 9.01E+00 2.77E+00 0 
Rosenbrock 
FSA 1.15E+01 1.52E+01 1.31E+01 6.88E-01 0 
(100 维 ) 
WPA 2.59E-11 5.26E-02 3.83E-02 8.48E-03 90% 
DWPA 3.23E-13 2.82E-10 2.39E-11 6.20E-11 100% 
PSO 5.04E+00 6.99E+00 6.08E+00 3.83E-01 0 
DE 1.92E+01 2.00E+01 1.98E+01 1.13E-01 0 
Ackley-1 
FSA 2.13E+01 2.15E+01 2.14E+01 3.42E-02 0 
(100 维 ) 
WPA 7.30E-13 6.36E-02 5.73E-02 1.76E-02 87% 
DWPA 8.88E-16 8.88E-16 8.88E-16 0 100% 
PSO 2.71E+00 3.29E+00 3.14E+00 1.42E-01 0 
DE 1.91E+01 1.98E+01 1.96E+01 1.13E-01 0 
Ackley-2 
FSA 2.13E+01 2.16E+01 2.14E+01 3.45E-02 0 
(200 维 ) 
WPA 1.52E-11 7.17E-01 2.39E-01 4.68E-01 75% 
DWPA 8.88E-16 8.88E-16 8.88E-16 0 100% 
人 对 于 60 维 的 多 峰 、 可 分 函数 Rastrigin，DWPA 和 WPA 10 
一 一 WPA 一 一 DWPA 
寻 优 成 功 ，DWPA 寻 优 精 度 高 于 WPA， 如 图 2 所 示 。 
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图 2 DWPA 与 WPA 的 Rastrigin 函数 平均 迭代 
g) 对 于 高 维 、 多 峰 的 不 可 分 函数 Griewank 逢 


较 稳定 , WPA 存在 陷入 局 


情况 而 使 得 寻 优 能 力 不 稳 定 ,如 图 3、4 所 示 。 


2000 
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部 最 优 的 


h) 对 于 高 维 、 


多 峰 的 可 分 函 


图 3 DWPA 与 WPA 的 Griewank 函数 平均 迭代 曲线 


数 Quadric 和 单 峰 的 可 分 函数 


Sumsquares-2，DWPA 和 WPA 寻 优 成 功 ， 且 DWPA 寻 优 精度 


高 于 WPA， 如 图 5 及 图 6 所 示 所 示 。 


ij) 对 于 高 维 的 多 
较 稳 定 ，WPA 存在 陷入 


峰 、 不 可 分 函 


数 Ackley，DWPA 寻 优 性 能 


局 部 最 优 的 情况 ， 且 随 着 Ackley 函数 
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维 数 的 增加 ， 其 陷入 局 部 最 优 的 可 能 性 增加 而 导致 寻 优 能 力 继 
续 下 降 ， 如 图 7、8 所 示 。 
ee 
函数 都 具有 较 好 的 寻 优 计算 能 力 , 与 WPA 相 比 , DWPA 对 
Re ents 数 的 求解 不 易 陷入 局 部 最 优 值 ， 算 法 鲁 棒 性 
也 得 到 加 强 。 
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图 6 DWPA 与 WPA 的 Sumsquares-2 函数 平均 迭代 曲线 
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结束 语 
为 了 提高 狼 群 算法 的 寻 优 性 能 ， 本 文 提出 一 种 基于 差分 进 


化 的 改进 狼 群 算法 (CDWPA), 在 游 走行 为 中 引入 探 狼 搜索 因子 ， 


使 得 算法 寻 优 能 遍历 整个 空间 ; 
过 程 ， 
条 时 


在 


在 奔袭 中 除去 人 工 狼 距离 计算 


引入 狐 狼 最 大 奔袭 次 数 和 头 狼 蔡 换 作为 两 个 行为 终止 
， 防 止 算法 陷入 无 限 循环 并 提高 计算 效率 和 算法 鲁 棒 性 ; 


所 攻 行为 中 引入 自 适 应 围攻 步 长 ， 使 得 距离 头 狼 相对 较 远 的 


人 工 狼 以 相对 较 大 的 步 长 向 头 狼 靠 近 
工 狼 以 较 小 步 长 对 头 狼 附近 


， 上 距离 头 狼 相 对 较 近 的 人 
区 域 进 行 更 加 精细 的 搜索 ， 兼 顾 算 


法 的 局 部 搜索 和 全 局 搜索 能 力 ， 并 在 改进 WPA 算法 的 基础 上 


加 入 差分 进化 的 变异 、 交 又 和 选择 
励 ， 提 高 算法 全 局 搜索 能 力 。 通 过 
5 种 算法 进行 了 性 能 对 比 测试 。 测 


操作 ， 使 算法 能 跳出 局 部 最 
13 个 测试 函数 对 DWPA 等 
试 结果 表明 ，DWPA 对 各 类 


压 


单 峰 /多 峰 、 低 维 / 高 维 、 可 分 /不 可 分 函数 都 具有 较 好 的 寻 优 计 
算 能 力 ， 其 在 求解 多 峰 、 高 维 、 不 可 分 函数 时 能 有 效 跳出 局 部 


最 优 ， 为 高 维 、 
和 高 效 的 解决 六 
如 何 对 游 走 、 
优 效果 需要 进 


等 复杂 函数 的 求解 提供 了 一 种 新 的 思路 
方 本 文 只 在 围攻 行为 中 引入 了 自 适应 步 长 ， 
奔袭 行为 中 的 固定 步 长 进行 改进 以 获得 更 佳 的 寻 
步 研究 。 
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